Excitation of Vortices in Semiconductor Microcavities 
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We predict that the transverse electric-magnetic polarization splitting of exciton-polaritons allows 
a simple technique for the generation of polariton vortices of winding number 2 in semiconductor 
microcavities. The vortices can be excited by circularly polarized light having no orbital angular 
momentum and observed as phase vortices in the opposite circular polarization or directly as linear 
polarization vortices. The prediction is explained by a simplified analytical linear model and shown 
fully with a numerical model based on the Gross-Pitaevskii equations, which includes the non-linear 
effects of polariton-polariton interactions. 
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Introduction. The creation and evolution of vor- 
tices in atomic Bose-Einstein condensates (BECs) has at- 
tracted much attention p], Q since the experimental real- 
ization of a dilute atomic BEC [3] just over a decade ago. 
It is no coincidence that vortices are at the heart of our 
understanding of superfluidity, forming spontaneously in 
type II superconductors [4] and superfluids 0,0). In 
atomic BECs vortices can be optically created using 
light modes with orbital angular momentum (Laguerre- 
Gaussian modes) and exploiting the recent technology of 
ultraslow light pulses to couple light and matter fields [7[ . 
Vortices written by light fields are stored in the atomic 
condensate and can later be re-written onto light fields. 

In semiconductor systems, the high temperature BEC 
of the half-light half-matter quasiparticles called exciton- 
polaritons has been reported recently [8[ . The condensa- 
tion takes place in a system of several quantum wells 
sandwiched between two distributed Bragg reflectors, 
commonly known as a microcavity. Many unique prop- 
erties of polaritons [9] arise from their non-parabolic in- 
plane dispersion (center left of Fig. [1]). No vortices of 
polaritons have been observed so far, although super- 
fluidity in polariton systems has been largely discussed 
theoretically 0, El • Here we show that a unique fine 
structure of polaritons in microcavities allows for a very 
simple mechanism of exciting vortices using conventional 
circularly polarized light beams having zero orbital an- 
gular momentum. 

It is well known that transverse electric (TE) and mag- 
netic (TM) normal modes of a cavity are non-degenerate 
for finite in-plane wave vectors. Polaritons also possess 
a polarization dependent TE-TM energy splitting [12[, 
which was previously proposed to allow the optical spin 
Hall effect [13| . We will show that this splitting can cause 
an initially circularly polarized polariton distribution to 
form a phase vortex in the cross-circular polarization, 
characterized by a phase singularity around which the 
phase changes by 47r. A vortex structure also appears in 
the linear polarization pattern. 

This linear effect is a general property of all planar 
cavities, not only those with quantum wells, and could 



be useful for creating light modes with angular momen- 
tum. Such light modes can have interesting applications 
to quantum information Q. Our next step is to take 
into account the many-body polariton-polariton interac- 
tions, which play a crucial role in a number of micro- 
cavity processes where large polariton populations are 
involved 0, [3. 

Polariton-polariton interactions can be accounted for 
using the zero-range interaction and mean-field approxi- 
mations that lead to the Gross-Pitaevskii equations. In 
Ref. [10[ the polarization independent spatial dynamics of 
interacting polariton systems was derived by numerically 
solving the single component (scalar) Gross-Pitaevskii 
equations. For us, the inclusion of polarization is not 
only important due to effects caused by TE-TM splitting; 
polariton-polariton interactions themselves were previ- 
ously shown to be spin-anisotropic, which led to effects 
such as self-induced Larmor precession of elliptical polar- 
izations and inversion of linear polarizations [17[ • 

Taking into account polarization, polaritons represent 
a two component Bose gas. Vortices in two-component 
bosonic systems are closely related to Skyrmions (l8j , 
which makes microcavities a potentially suitable solid 
state system for Skyrmion observation. 

Linear Model. We will first outline a simple analytic 
model to show how vortices can be generated in micro- 
cavities. First we consider a lower branch polariton state 
characterized by the in-plane wavevector k rotated by an 
angle <p with respect to the x— axis. A circularly polarized 
state is a linear combination of TM and TE eigenstates, 
given by: 
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The evolution of the initial state in the circular basis is: 
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where uutm and uote are the frequencies of the two eigen- 
states and r accounts for the life-time of polaritons. The 
k-dependent energy splitting between uutm and uote is 
shown in the center right plot of Fig. [TJ From Eq. (|3|) 
we see that the phase of the cross-circularly polarized 
component (-0_) is angle dependent. So far we have only 
considered a single state with wavevector k. To obtain 
the field distribution from multiple wavevector states we 
should take the convolution integral with a pump distri- 
bution, /(fc, 0, i): 
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where # is the angle between the radius vector r and the 
x— axis. For simplicity we consider excitation by a delta 
function pulse f(k,<p : t) = AS(k — k p )S(t). This gives: 
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where J n denotes the nth order Bessel function of the 
first kind. The amplitude and phase of and are 
plotted in Fig. [TJ It shows that the phase of the cross- 
circular component changes by 4tt around the origin. In 
other words we have found a vortex of winding number 2. 
Note also the role of oscillatory terms in Eqs. ([3]) and ([5]): 
they describe the beats between the two spin states of 
polaritons. Initially, all polaritons lie in the spin- up state 
and there is no vortex. After a quarter of the period T = 
— all polaritons lie in the spin-down state and the 
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vortex appears. Another quarter period later polaritons 
return to the spin- up state and the vortex disappears. At 
the intermediate time t = ^ the polaritons are linearly 
polarized, forming a polarization vortex as sketched in 
the bottom left plot of Fig. [TJ 

Non-Linear Model. In this section we use a numeri- 
cal model to include the effects of polariton-polariton in- 
teractions and the finite extent of any real pulse in space 
and time. We will use the Gross-Pitaevskii equations, 





FIG. 1: {colour) (Top:) Plots of the absolute values of ^+ 
(left) and \£_ (right) for k p = 100mm -1 . Note that the gray- 
scales are not the same; the cross-circular component is much 
weaker (depending strongly on the value of the TE-TM split- 
ting at k p ). (Center Left:) The dispersion of photons and 
excitons (gray) and polaritons (black) for an exciton-photon 
coupling energy (vacuum field Rabi splitting) of 5.1meV. The 
TE-TM splitting of the photon is given by the different effec- 
tive masses of TE and TM branches. We chose rat = 10 -5 ra e 
and mi = 0.95rat where ra e is the free electron mass. (Center 
Right:) The parabolic TE-TM splitting of the photon causes 
an TE-TM splitting of the lower (solid) and upper branch 
(dashed) polaritons. (Bottom Left:) The linear polarization 
vortex. Arrows indicate the linear polarization along the cir- 
cle k — k p after a time t — 2 ^ t ^_ Ute ^ . (Bottom Right:) 
The phase of \£_ changes by 4tv around the origin. 



which are a familiar tool for the treatment of dilute Bose- 
gases [l9|, [lo| . For a coherent ensemble of polaritons the 
set of coupled Gross-Pitaevskii equations is written: 
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(we sum over the index j). Equations (|6l[7j) are obtained 
by the mean field treatment, valid at low temperatures 
and densities where the thermal and quantum depletion 
of the condensate has only a small effect. Xi an d <$>% are 
wavefunctions representing the excitonic and photonic 
components of the upper and lower branch polar itons. 
The indices z,j can take the values + or — , represent- 
ing the two circular polarizations. The operators T x and 
are the kinetic energy operators of uncoupled exci- 
tons and photons, which have parabolic dispersion with 
effective masses m x and m§\ 
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We have chosen an energy scale such that the energy of 
the photonic component is zero at zero in-plane wavevec- 
tor and consider the case where there is no detuning be- 
tween the exciton and photon modes. The parameter 
£ allows for a (parabolic) TE-TM splitting of the pho- 
tonic component in the microcavity (we neglect the much 
smaller splitting of the excitonic component). 

The coupling between excitons and photons is de- 
scribed by the Rabi splitting of the cavity, 2V, which 
is related to the quality factor and the exciton decay 
rate @, [2l[. The non-linear terms in Eq. (|6]) describe 
Coulomb interactions between the excitonic components 
where a\ describes the parallel spin configuration and (\2 
describes the anti-parallel configuration. In Eq. (0) the 
term fa (x, i) represents the optical pumping of the micro- 
cavity. For a coherent Gaussian pump pulse exciting the 
lower polariton branch at the wavevector ko and energy 
E p at time t p the pump in reciprocal space is: 
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Ai represent the amplitudes of the two circular polariza- 
tions of the pump pulse. Ak and At define the width 
of the pump in reciprocal space and time. The frac- 
tion T/ (E LP (k) - E p - irf) is the linear susceptibility 
of a system of Lorentz oscillators where T is the homoge- 
neous oscillator (HWHM) linewidth [22[. In Eq. (7J) the 
pumping term t) is the inverse Fourier transform of 
/(fc, t). ELp(k) is the bare dispersion of the lower polari- 
ton branch, which is calculated by finding the eigenval- 
ues of the linear Hamiltonian (without interaction terms) 
in the exciton-photon basis [23[. Finally, r x and are 
the lifetimes of the excitonic and photonic components, 
which account for the inelastic scattering and radiative 
decay of polaritons respectively. 
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FIG. 2: (colour) Plots of the absolute values of the pho- 
tonic wavefunction squared, |0;| 2 , at different times for the 
co-circular (left) and cross-circularly polarized components 
(right). The pulse was centered at t p = Ops and had widths: 
At = lOps and AL = 2/Ak = 10/mi. Other parameters: 
E p -4.75meV, V 5.1meV, T O.lmeV, ( 0.025, 
r v — 50ps, — 3.3ps, ai — 5 x 10 _5 meV mm 2 , — — O.lai, 



'x 

771 



: 0.22m e and — 10 m e 




3.14 



-3.14 




3.14 



■3.14 



FIG. 3: (colour) Plots of the phase of the photonic wavefunc- 
tion, arg(0i), at 7.5ps after the arrival of the pulse for the 
co-circular (left) and cross-circular polarizations. 



Equations (|6]) and Q completely determine the dy- 
namics of interacting polaritons once initial wavefunc- 
tions and pumping terms are defined. We acknowledge 
that the polariton dynamics from a linear pulse pump 
was recently studied [11] using similar equations (and 
our model was able to reproduce these earlier results). In 
this paper we focus on the case of excitation with a circu- 
larly polarized pump pulse (A+ = 1, A- = 0) with broad 
spread in wavevector tuned slightly (0.35meV) above the 
bottom of the lower polariton branch. Although the 
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pump is centered at ko = (normal incidence) it is 
resonant with the lower polariton branch along a ring 
(k = 435mm -1 ) in reciprocal space with a radius set by 
the pump energy. We numerically solved Eqs. (|6]) and 
(|7j) to find the results plotted in Figs. [2]and[3l We show 
the behaviour of the photonic wavefunctions, which are 
directly accessible experimentally. The excitonic wave- 
functions demonstrate an identical behaviour. 

The appearance of expanding rings in the distributions 
should be expected since we are exciting a ring in recipro- 
cal space. In the cross-circular polarization there is again 
a phase singularity and vortex of winding number 2 (in 
agreement with the analytical linear model). 



coherent polariton scattering, introduced via dissipative 
terms in Eqs. ([7j) and (|6]), does not affect the appearance 
of the vortices, whilst it does influence the intensity of 
polarized light emitted by the cavity. 

Conclusion. We calculated the dynamics of an in- 
teracting polariton condensate, which was excited by a 
circularly polarized Gaussian pump pulse tuned above 
the bottom of the lower polariton branch. The presence 
of TE-TM splitting in the microcavity allows coupling 
into the cross-circular polarization in which the phase 
shows the profile of a vortex of winding number 2. The 
appearance of a linear polarization vortex and spreading 
ring-like structures is also derived. 
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FIG. 4: (colour) The intensity of the photonic wavefunction, 
|0i| 2 , in the co-circular (left) and cross-circular (right) polar- 
izations for different values of a\ at 15ps after the peak of the 
excitation pulse. We keep the relationship ct<2 — — O.lai. 

The expanding polariton rings are influenced signifi- 
cantly by the polariton-polariton interactions as shown 
in Fig. [4] where we plot cross sections of the wavefunc- 
tions along the x-axis for different values of the scatter- 
ing parameter, ol\. In the absence of polariton-polariton 
interactions (a\ — 0) a ring is observed in both circu- 
lar polarizations. As ol\ increases the intensity of the 
wavefunctions increases slightly (formally an increase in 
the interaction constants is identical to an increase in 
the pump amplitude) but as the interaction constants 
increase further the intensity decreases, and the wave- 
functions spread due to the repulsion of polaritons in 
real space. Also note that as ol\ increases we observe a 
central peak in the co-circular distribution instead of a 
ring. This is due to the larger (time-dependent) blueshift 
of the system when ol\ is increased, which results in the 
resonant excitation of a smaller ring in reciprocal space. 
In real space we also observe the appearance of darker 
rings within the brighter ring in cross-circular polariza- 
tion as ol\ is increased from zero. Such structures seem 
similar to dark ring solitons, which are known to appear 
in the equilibrium scalar Gross-Pitaevskii equations and 
similar non-linear Schrodinger equations [24Ll25|]. The in- 
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